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(54) Probe sequences intercalating method with information data for the estimation of 
communication channel characteristics and correct reception of data sequences 



(57) A method for estimation of a transmission 
channel and tor transmission and reception of digital 
signals over the channel comprises the steps of trans- 
mission of data sequences intercalated with probe 
sequences, calculation in reception of a channel char- 
acteristics estimation for probe sequences received and 
interpolation of the channel estimation for the data 
sequences. After obtaining the interpolated estimation 
the transmission and reception method comprises the 
additional phase of recognition of received information 
data by equalizing the data sequences with the interpo- 



lation result. An apparatus employing the method has at 
one end of the channel a device (12) for intercalation of 
probe sequences between data sequences and at the 
other channel end a discriminator (16) which separates 
the probed and data sequences. An estimator (17) esti- 
mates the channel for the probe sequences and an 
interpolator (18) produces the desired interpolated esti- 
mation. An equalizer (19) equalizes the data sequence 
on the basis of the interpolated estimation. 
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Description 

[0001] The preserrt invention relates to a method and system in accordance with this method tor the implementation 
of numerical connections in particular through channels dispersive in frequency and rapidly varying in time in addition 
5 to being noisy. 

[0002] The employment of radioetectrical channels for high-speed numerical connections is limited by the presence 
of intersymbolic interference (ISI) and propagation fading with characteristics generally varying in time and caused by 
the multiplicity of paths present between the transmitter and the receiver. In this situation the signal received is 
expressed in the equivalent model (complex) in time-discrete base band (symbol Ts time sampled), as: 

10 

L-1 

y(0= X g(i;m)a(i-m)+v(i) 

m-0 



15 

where a(i) is the numerical sequence transmitted (e.g. PSK, QAM, OMSK equivalent), v(i) is the added noise (thermal 

noise plus external interference) and g(i;m). m=0 L-1 are the L coefficients of the time-variant impulsive response of 

the transmtssive channel including the transmission and reception filters. 

[0003] In particular the temporal variability of the transmission channel can be induced by the mobility of the transmit- 

20 ter and/or receiver as happens in mobile radio connections where it increases with the speed of the vehicle or by the 
fluctuations in time of the altitudes of the strata of the atmosphere as is found in HF radio bridges. This variability can 
be expressed by the coherence time tc which is defined as the interval of time within which the channel does not vary 
significantly from a statistical view point or, equivalent!* in terms of Doppler spread Bd=1Ac. Normalizing tc with 
respect to the symbol time Ts gives the number Lc=tc/Ts of consecutive symbols for which the channel can be conskJ- 

25 ered nearly constant in time. In HF or mobile radio applications typical values of the product of Bd-Ts go from 10 4 
(slowly varying channel) to 10~ 3 (tast channel) up to a maximum of 5x1 0" 3 or 10" 2 (e.g. for high-speed vehicles like fast 
trams or planes); contrariwise, for Bd - Ts less than 10" 4 the channel can be considered virtually static. 
[0004] To obtain satisfactory performance in terms of Bit Error Rate (BER) the presence of ISI makes necessary the 
use of non-linear equalization techniques. Among these the Maximum Likelihood Sequence Estimation (MLSE) algo- 

90 rrthrn displays reduced computational complexity thanks, for example, to its implementation by the Viterbi algorithm 
(MLSE-VA) and lor virtually static channels it provides broadly optimal performance since the decision on a certain sym- 
bol received takes place on the basis of the symbols either preceding or following it. The MLSE operates on the basis 
c* a decision delay of 0 samples received where D is generally chosen several times the length L of the impulsive 
response of the channel. The efficiency of the MLSE solution in its different versions known in the technical literature is 

35 such that it is commonly employed to equalize frequency-selective channels. e.g. in GSM receivers. 

[0005] When the transmission channel varies in time the MLSE must be assisted by a channel estimator which at the 
beginning of the transmission makes the channel estimation and then during reception of the data continuously updates 
the estimation which is thus supplied continuously to the MLSE equalizer. This solution is termed adaptive MLSE or 
AMLSE. 

40 [0006] The channel estimator typically employs the data decisions at the MLSE output and is commonly implemented 
by an algorithm of the gradient or a Kalman filter (see for example sec. 6.8 of the book by J.G.ProaHs, Digital Commu- 
nication, second edition, McGraw-Hill, 1 989) which is the optimal solution. But in any case the updating takes place with 
a delay equal to the decision delay D. 

[0007] rt is known that the AMLSE equalizer loses much of its efficiency when the channel varies rapidly in time 
<5 because the channel estimation is constantly delayed by D symbols with respect to the data and in this interval the 
channel may have varied in a not negligible manner. To obviate this problem it is posstole to supply the channel estima- 
tor with low^delay Trial* decisions to then make the final decisions with a greater decision delay. This mitigates but 
doesn't solve the above mentioned problem because the trial decisions are less reliable than the final ones. Another 
possible solution is the one based on the so-called Per-Survivor Processing (PSP) which employs an MLSE equalizer 
so in which a different channel estimation is performed for every possible surviving path. This gives significant improve- 
ments in performance at the cost however of considerable computational complexity caused by the high number of 
channel estimators needed. 

[0008] The general purpose of the present invention is to obviate the above mentioned drawbacks by making available 
a method and a transmission system applying this method for transmission and reception of numerical signals which 
55 would allow high efficiency even in the case of noisy, frequency-dispersive channels rapidly variable in time. 

[0009] In view of this purpose H was sought to provide in accordance with the present invention a transmission and 
reception method tor digital signals comprising information data propagated in a communication channel in which are 
present the following steps. Implement a frame comprising information data sequences to be transmitted spaced by 
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probe sequences consisting of known data, transmit the frame in the channel, receive the frame from the channel, cal- 
culate a channel characteristics estimation for probe sequences contained in the received frame, perform for data 
sequences an interpolation of the channel estimation found for the probe sequence, and perform acknowledgement of 
information data received equalizing the data sequences by the use as a channel characteristics estimation for the data 

5 sequences of the result of the interpolation. 

[0010] Again in accordance with the present invention it was sought to provide a method for estimation of the charac- 
teristics of a communication channel during transmission over it of information data sequences, comprising the follow- 
ing steps. Intercalate with the transmitted data sequences probe sequences consisting of known data, calculate the 
channel characteristics estimation for probe sequences, interpolate the calculated estimation, and employ the interpp- 

w lation result as a channel characteristics estimation during reception of information data sequences. 

[001 1] It was also sought to provide an apparatus for transmission and reception of digital signals comprising infor- 
mation data propagated in a communication channel characterized in that it comprises a transmitting apparatus and a 
receiving apparatus at the two ends of the channel with the transmitting apparatus comprising in turn an intercalating 
device receiving at input the information data to be transmitted and outputing a frame made up of information data 

is sequences to be transmitted intercalated with probe sequences consisting of known data and a transmitter for sending 
the frame over the channel with the receiving apparatus comprising in turn a receiver for reception of the frame from the 
channel with a discriminator separating the probe sequences and the data sequences from the received frame, an esti- 
mator receiving at input the probe sequences and supplying at output the channel estimation characteristics for the 
received probe sequences, an interpolator receiving at input the estimation output by the estimator and producing by 

20 interpolation the desired channel characteristics estimation for the data sequences and an equalizer receiving the esti- 
mation output by the interpolator and equalizing on the basis thereof the data sequence extracted from the discriminator 
to output the correct sequence of transmitted data. 

[0012] To clarify the explanation of the innovative principles of the present invention and its advantages compared 
with the prior art there is described below with the aid of the annexed drawings a possible embodiment thereof by way 
25 of non-limiting example applying said principles, in the drawings: 

fig. 1 shows transmitted frame structures provided in accordance with the present invention. 

figs. 2 and 3 show graphs of the channel estimation produced in accordance with two possible variants of the 

present invention, and 

30 - fig. 4 shows a block diagram of a possfcle digital signal transmission and reception system applying the present 
invention. 

[001 3] In accordance with the present invention with the information data or data sequences transmitted are interca- 
lated probe sequences or probe data as shown by way of example in fig. 1 . The probe data are known a priori and can 
35 be intercalated in a variety of manners depending also on the structure of the information data which it is desired to 
transmit on the transmission channel coherence time and on the reception noise level. 

[0014] For example, in fig. 1a transmission of the data consisting of PSK or QAM symbols is structured like a contin- 
uous sequence of symbols in which the sequences of information data D consisting of a number of symbols defined by 
Ldata are alternated with probe sequences S consisting of a number of symbols defined by Lprobe. Each probe 

40 sequence acts either as a preamble for the data following it or postambie for the data preceding it. A data sequence and 
the probe sequence preceding it constitute a "frame" formed of a number of Lfra symools equal to Lfra=Lprobe+Ldata . 
The durations in time of the frame, probes and data between the probes are thus respectively equal to Tfra=Ts x Lfra , 
Tprobe=Ts x Lprobe and Tdata=Tfra-Tprobe with Ts= time symbol. A block of symbols can be seen as the assembly 
of a block of data and the immediately preceding and following probe sequences. 

45 [001 5] Successive blocks of data thus defined are therefore superimposed and have in common the probe sequence 
between them. 

[0016] in fig. lb transmission takes place in independent blocks consisting each of a preamble of probe data Spre, 
information data D and a postambie of probe data Spos. The number of preamble symbols Lpre and postambie sym- 
bols Lpost can be equal mutually and to Lprobe. 
so [001 7] In fig. 1 c the transmission takes place in consecutive blocks, i.e. one or more consecutive blocks are transmit- 
ted separated by intermediate probe sequences (termed "midambles") Smid to form independent "packages" starting 
with a preambular probe sequence and a postambular probe sequence. The number of midamhular symbols can be 
equal to the number or preambular and postambular symbols. 

[0018] For the sake of simplicity, by the words "preambular and postambular probe sequences" is meant the pair of 
55 probe sequences between which there is an information data sequence regardless of the fact that they are in reality 
intermediate probe sequences or probe sequences which are simultaneously preambular for the following data 
sequence and postambular for the preceding data sequence. 

[0019] In accordance with the present invention, upon reception of the signal the channel is estimated for only the 
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probe sequences. The estimation thus found is interpolated to obtain a channel estimation for the data sequences. The 
interpolation can be advantageously linear. 

[0020] The channel estimation thus found is employed to perform equalization of the information data output by the 
receiver. 

5 [0021] A transmission system in accordance with the present invention ts shown in the block cfiagram of fig. 4 and 
indicated generically by reference number 10. The transmitting part 1 1 receives a numerical sequence of information 
data Din to be transmitted (e.g. PSK. QAM, GMSK-equivalent etc.) . An intercalated 1 2 introduces into the Din sequence 
the probe sequences S to provide a predetermined frame structure to be transmitted (continuous, with independent 
blocks, consecutive blocks etc). The output of the intercalates is started at a signal transmission block 1 3. Intercaiator 

ro and transmission block are known in themselves and reac% imaginable for those skilled in the art and are according 
not further described or shown. 

[0022] Trie receiving part 1 4 comprises in turn a receiving block 15 homologous with the transmitting block 1 3. The 
received signal is sent to a discriminator 16 which separates the probe and data sequences. The probe sequences are 
sent to an estimator 17 which outputs the channel estimation for the received probe sequences. The estimation is sent 
is to an interpolator 18 which produces the desired channel estimation for the data sequences. The latter estimation is 
sent to an equalizer 1 9 which on the basis thereof equalizes the received data sequence and outputs the correct data 
Dout to the transmission system. 

[0023] It has been found advantageous that the received sequence be encoded with channel coding of the known 
Trellis type or the convolutional type. It has also been found advantageous that the received sequence be subjected in 
20 transmission to codification of the differential type or with interlacing. 

[0024] The algorithm for channel estimation during the probe sequence should be chosen from among the estimation 
algorithms with the characteristic of converging rapidly on the estimated value in addition to giving rise to a small esti- 
mation error of course. In this manner the Lprobe length is kept small and consequently high Ldata/Lfra transmission 
efficiency is achieved. 

25 [0025] For this purpose any fast estimator, e.g. of the Kahman type, can be used. But the Kalman-type non-linear 
recursive adaptive estimator described in European patent application n° 97 202 349.3 filed 25 July 1997 in the version 
related to known input data has been verified as being particularly advantageous. With this estimator the best perform- 
ance was achieved in terms of estimation precision and convergence speed. 

[0026] Briefly, the estimator is based on the following recursive procedure which is valid in general for complex mod- 
30 ulations of the QPSK or GMSK-equivalent type. Trie steps of the procedure are as follows. 

0) Let x(i)=[a(i) a(i-L+1)] T be the L-dimensional vector (complex) of the "status- of the channel at the Hh instant 

known to the receiver because the transmitted data a(i) belong to the probe sequence and let g(i)»[g(i;0) g(i;L)] T 

be the vector of the channel coefficients at ith instant. Here and hereinafter v T indicates transposition ol vector v. 

35 1) in step ith calculate as g A (i/i-l)=A g^-ITM) the prediction (one step) of g(i), i.e. the estimation of g® based on 
the observations up to the step (M)th, where g*(i-1/i-1) is the "filtered" estimation of g(i) found in the preceding 
step (i-l)th while A is the channel status transition matrix, real and having dimensions LxL, whose elements are 
found by modeling the channel process as a first order autoregressive process (vectorial, with dimensions L). The 
"filtered" estimation is initialized in the initial step (i.e. for the beginning of the probe sequence) as described below. 

40 2) Calculate recursively the LxL matrix of the covariances of the channel prediction error as: 

Sg(i/i-1)=A[l UL -K(i-1)x T (i-1)]Sg(i-1/i-2)A T + +2[R fl (0)-AR g (0)A T ] , where Rg(0) is the covariance matrix of the 
channel process g(i), K(i-1 ) is the f flter gain found in the preceding step and l^ is the identity matrix with dimen- 
sions LxL 

3) Calculate the filter gain as: 

45 

K(i)=[Sg(i/i-1)M x(i)/{N 0 +x(i) T M H Sg *(i/i-1)Mx(n)} 

where M is the "mapping" Lxisj matrix associated with the N possible statuses of the channel while No is the noise 
50 power spectral density. Here and below the notation B* indicates the conjugate of a matrix or vector B. 

4) Calculate the "filtered" estimation of g(i) as: 

g^W-^-iHKOJtyO^^-i) T MxC0] 

55 

5) If the end of the probe sequence is reached g*(W) constitutes the channel estimation associated therewith; oth- 
erwise increase the index i and return to step 1). 
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[0027] It has been found advantageous to transmit the preambles in the simple BPSK format, i.e. on two levels, what- 
ever be the modulation format (e g. QPSK) employed for transmitting the information data. Indeed, in this manner the 
channel estimation is just as good as with QPSK which however is more onerous from a computational viewpoint. 
[0028] Bigger modulation formats give rise to less accurate channel estimations. 
5 [0029] In addition, the precision in the channel estimation at the end of the probe sequence can be improved by 
repeating the above described estimation several times always employing the same probe sequence but starting each 
time from the estimation found at the end of the preceding iteration. In this manner there can be obtained convergence 
of the channel estimation algorithm even when Lprobe is small and a single iteration is insufficient for achieving conver- 
gence. 

io [0030] As mentioned above, the information data channel estimation is found as the probe sequence channel estima- 
tion interpolation. 

[0031] To obtain this, the channel is estimated on the basis both of the preceding probe sequence and of the probe 
sequence which follows an information data sequence. 

[0032] The behavior of the channel for the Ldata information symbols present between~the two probes is recdn- 
15 struct ed by means of an interpolation algorithm. It has been found advantageous to employ a simple linear interpolation 

algorithm as illustrated graphically in fig. 2 in which is shown a possible trajectory (drawn in solid lines) of the real part 

Re{g(i;m)} or the imaginary part lm{g(i;m)} of the generic channel coefficient g(i;m). For the sake of simplicity this real 

or imaginary part is still indicated here with g(i;m) without thereby introducing ambiguity or confusion. 

[0033] The channel estimator follows the trajectory for the preambular and postambular probe sequences (in the fig- 
20 ure possible estimation trajectories are drawn in broken lines) and gives at the end thereof the estimations g*pre and 

g A post respectively. With reference to fig. 2, using a linear interpolation the trajectory of g(i;m) (in broken lines in the 

figure) is reconstructed as: 

- v (g"post(i-i^e)^^pre[Ldata^Lpre-(M_pre)]) 
2S *V'™> (Ldata+Lpre) 



where i __pre is the preamble end instant and therefore the postamble end instant is Lpre+Ldata+Lpost . 

[0034] Since the above described linear interpolation tends to lose precision far from the g*pre and g A post estima- 

30 tions, i.e. in the center of the data block, it is possible to repeat once more the channel estimation algorithm starting 
however from the g"post estimation and travelling over the postamble in the opposite direction to its initial point. The 
estimation thus found g^st 1 is employed for the linear interpolation instead of g A post in a manner similar to that 
described above. This is shown graphically in fig. 3. Once the channel estimation has been found by interpolation for 
the data sequence recognition of the data can take place by an appropriate version of the known algorithm MLSE-VA 

35 with decision delay D by using said channel interpolated estimation symbol by symbol. In particular the metrics for the 
VA decision trellis stage for the ith instant are calculated step-by-step using the channel estimation for that same 
instant. The VA departs from a known initial status calculated on the basis of the last known preamble symbols and ter- 
minates with a known final status calculated on the basis of the first known postambular symbols. 
[0035] Channel estimation for the probe sequences, interpolation of the channel estimation for the information data 

40 and recognition of the data are repeated iteratively for reception of the following blocks. If transmission is in continuous 
sequence as shown in fig. la or with consecutive blocks separated by midambles as in fig. 1c the g A post estimation 
constitutes the g A pre of the following block which therefore does not need to be recalculated. With every iteration the 
channel is therefore estimated for the postamble while the preamble estimation is already available. This helps 
decrease the computational complexity of the method in accordance with the present invention. 

45 [0036] At the beginning of the probe sequence it is necessary to initialize the channel estimation and the correspond- 
ing error matrix Sg(iA-l). In the case of continuous sequence transmission as shown in fig. 1 a or in consecutive blocks 
as in fig 1c and if Tdata is not greater than the channel coherence time tc the recursive channel estimator can advan- 
tageously be initialized with the channel estimation found for the preceding probe sequence which constitutes a good 
initial approximation and gives rise to a fast convergence (as may be seen in fig. 2) so that Lprobe can be chosen rela- 

so tively small. 

[0037] For the first probe sequence transmitted the channel estimator is initialized at zero and it can be expected to 
employ for this probe sequence a probe sequence a bit longer than the following probes. 

[0038] Again in the case of continuous sequence or consecutive block transmissions and for Ldata so high it is 
assumed that the channel varies considerable from one probe to the next it was found more appropriate to initialize the 
55 recursive channel estimator at zero at the beginning of each probe. The probe sequences must then have a length 
Lprobe sufficient to ensure convergence possibly even using the expedient of repeating the estimation several times 
always employing the same probe sequence but starting each time from the estimation found at the end of the preced- 
ing iteration. 
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[0039] in the case of a transmission with independent data blocks like that shown in fig. 1b the channel is in general 
much different from one block to another among other reasons because different blocks could be transmitted on bands 
with different frequencies as in so-called frequency hopping, and the channel estimator must be initialized at zero for 
the preamble. In this case it is possWe to provide the use of a preamble slightly longer than the postamble. 

s [0040] In all the cases described above the co variance matrix of the prediction error Sg(i/i-1 ) is initialized as equal to 
the covariance matrix of the error between the initial value assumed for the channel and the true value of the channel. 
[0041 ] However the structure of the frame transmitted is composed - formed of probe sequences and information data 
- it has been found particularly advantageous that Lprobe be rather small compared with Ldata so that the transmission 
efficiency Ldata/Lfra is between 30% and 90% or even more and in particular al least greater than 50%. 

io [0042] But Lprobe must be high enough to alio* convergence of the iterative estimation algorithm described above, 
possibly repeating the estimation several times always using the same probe sequence but starting each time from the 
estimation found at the end of the preceding iteration, and to supply a good channel estimation when there is noise. On 
this subject it should be clarified that since the channel estimator described above by way of example requires that the 
status of the channel be known, in reality the estimation can departs 

15 which a number of symbols equal to Lprobe-L+1 instead of Lprobe are thus used. 

[0043] Ldata must not be too big, otherwise the channel could vary considerably from one preamble to the next and 
thus reduce the precision of the following interpolation. 

[0044] Since the data decision takes place only after having received the postamble, i.e. with a delay of 
Tdata+Tprobe, to which must be added the decision processing time, Ldata must not be too high so as to keep this 

20 delay within the system specifications. 

[0045] At this point is its clear that the preset purposes have been achieved allowing transmission and reception with 
excellent results of digital signals which propagate in a communication channel with multiple time-variant paths which 
give rise to time and frequency distortions. The transmission strategy provided in accordance with the present invention 
is based on the observation that even for very high Bd • Ts values (e.g. up to 10" 2 ) the properties of the channel vary 

25 slowly with respect to the symbol frequency. Consequently it is possible to measure periodically on the basis of the 
probe sequences the values assumed by the channel coefficients g(i;m) and then reconstruct its behavior between one 
probe and the next (i.e. for information data transmission) by means of an appropriate interpolation. 
[0046] If Ldata is not high, the linear interpolation supplies a good channel estimation along the data sequence 
because, virtually, the curve which constitutes the channel trajectory is approximated by a succession of rather short 

30 rectilinear segments. 

[0047] In addition, it is efficient not to use the data decisions which in the presence of decision errors reduce the pre- 
cision of the estimation in the prior art systems. Not using the information data decisions there is also found consider- 
able reduction of the computational load with respect to the case of an adaptive equalizer performing the channel 
estimation along the entire received sequence, including the information data. 
35 [0048] The computational complexity of the channel estimator can be considered virtually the same as the classical 
Kalman's estimator, in its version for real data (but complex channel) rf the modulation is the BPSK type or, in any case, 
rf the BPSK format is employed, i.e. two level, for the probe sequences. 

[0049] The computational load is in reality less because the channel estimator is active only for the probe sequence. 
[0050] The complexity of the MLSE equalizer with delay D is the conventional complexity of the VA. The method and 
40 device in accordance with the present invention are therefore easy to implement using normally available circuits and 
technical solutions. 

[0051] Naturally the above description of an embodiment applying the innovative principles of the present invention 
is given by way of non-limiting example of said principles within the scope of the exclusive right claimed here. 
[0052] For example the estimator could be different from that described as the interpolation technique used could be 
45 different even though both were found advantageously useful in the above described form. 

Claims 

1 . Method for transmission and reception of digital signals comprising information data which propagate in a commu- 
so ni cation channel inducting the following steps: 

a) provide a frame comprising information data sequences to be transmitted spaced by probe sequences 
made up of known data; 

b) transmit the frame in the channel; 
ss c) receive the frame from the channel: 

d) calculate a channel characteristics estimation for the probe sequences contained in the received frame; 

e) for data sequences perform an interpolation of the channel estimation found in accordance with step d) 
above; and 
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f) perform recognition of information data received equalizing the data sequences by using the interpolation 
result as the data sequence channel characteristics estimation. 

2. Method in accordance with claim 1 characterized in that the frame comprises data blocks made up of an tnforma- 
5 tion data sequence preceded by a preambular probe sequence and followed by a postambular probe sequence. 

3. Method in accordance with claim 2 characterized in that each data block has postambular probe sequences con- 
stituting preambular probe sequences for an immediately following data block. 

io 4. Method in accordance with claim 3 characterized in that the data blocks are organized in block packages with each 
package being separated in time from the preceding and following packages. 

5. Method in accordance with daim 1 characterized in that the interpolation is a linear interpolation performed on the 
basis of the estimation of probe sequences immediately preceding and following the data sequence. 

15 

6. Method in accordance with daim 5 characterized in that the channel estimation for the probe sequence following 
the data sequence is found traversing said probe sequence in the opposite direction. 

7. Method in accordance with claim 5 characterized in that the trajectory gfi;m) of the generic channel coefficient for 
20 a data sequence is interpolated linearly as: 

f . ( g^pos^i^re^^pretLdata^Lp re-ri-i^re)]) 
91 ' ' (UJata+Lpre) 

where g*pre and g*post are the estimations calculated for the preambular and postambular probe sequences 
respectively before and after the data sequence. Lpre is the end of the preamble instant Ldata is the length of the 
data sequence, and Lpre is the length of the preamble. 

30 8. Method in accordance with claim 1 characterized in that the channel characteristics estimation for the probe 
sequences is performed by means of an adaptive non-linear recursive Kalman filter. 

9. Method in accordance with claim 8 characterized in that the Kalman filter is based on the recursive procedure of: 

35 a) In step rth calculate as g A (i/i-l )=A g A (i-1/i-1) the prediction (one step) of g(i), i.e. the estimation of g(i) based 

on the observations up to the step (i-1)th, where g A (i~1/M) is the "filtered" estimation of g(i) found in the pre- 
ceding step (i-1)th while A is the channel status transition matrix, real and having dimensions LxL, whose ele- 
ments are found by modeling the channel process like a first order autoregressive process (vectorial, with 
dimensions L); 

40 b) calculate recursively the LxL matrix of the covariances of the channel prediction error as: 

Sg(i/i-l)=AIl LxL -K(M)x T ^ Rg(0) is the covariance matrix of 

the channel process g(i), K(i-1 ) is the filter gain found in the preceding step and l LxL is the identity matrix with 
dimensions LxL; 
c) calculate the filter gain as: 



25 



45 



K(iMSg(i/M)lv1 *x(i)/{N 0 +x(i) T M H Sg *(i/i-1)Mx(n)} 



where M is the "mapping" LxN matrix associated with the N possfixe statuses of the channel while Nq is the 
so noise power spectral density; 

d) calculate the "filtered" estimation of g(i) as: 

g*(i/i)=g*(W-1 )+K(i)[y(i)-g A (i/i-1) T Mx(i)] 



55 



e) if the end of the probe sequence is reached, g^i/i) constitutes the channel estimation associated therewith; 
otherwise increase the index i and return to step a). 
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10. Method in accordance with daim 9 characterized in that if the duration in time of the data sequenoe Tdata is not 
greater than the channel coherence time tc the recursive procedure is initialized with the channel estimation found 
for a preceding probe sequence. 

5 11. Method in accordance with claim 9 characterized in that the recursive procedure is initialized at zero at the begin- 
ning of each probe sequence. 

12. Method in accordance with daim 1 characterized in that recognition of the information data is performed by means 
of algorithm MLSE-VA with decision delay D using said channel interpolated estimation symbol by symbol. 

10 

1 3. Method in accordance with claim 1 2 characterized in that the metrics for the decision trellis stage of the VA corre- 
sponding to the rth instant are calculated step-by-step using the channel estimation for that same instant with the 
VA departing from a known initial status calculated on the basis of the last symbols of a probe sequence immedi- 
ately preceding the data sequence and terminates with a known final status calculated on the basis of the first sym- 

15 bote of a probe sequence immediately following the data sequence. 

14. Method in accordance with claim 1 characterized in that the received frame is encoded with trellis channel coding. 

1 5. Method in accordance with daim 1 characterized in that the received frame is encoded with a convolutional channel 
20 coding. 

16. Method in accordance with claim 1 characterized in that the frame is subject in transmission to differential encod- 
ing. 

25 1 7. Method in accordance with claim 1 characterized in that the frame is subject in transmission to inter lacing. 

18. Method in accordance with daim 1 characterized in that the probe sequences are modulated with binary modula- 
tion and the data sequences are modulated with higher order modulation. 

30 19. Method in accordance with claim 1 characterized in that the channel estimation at the end of a probe sequence is 
found recursively by repealing several times the channel estimation tor the probe sequence and always employing 
the same probe sequence but departing each time from the estimation found at the end of the preceding iteration. 

20. Method for estimation of the characteristics of a communication channel during transmission over it of information 
35 data sequences comprising the following steps: 

intercalate probe sequences consisting of known data with the data sequences transmitted; 

- calculate the channel characteristics estimation for probe sequences; 
interpolate the calculated estimation; and 

40 - employ the result of the interpolation as a channel characteristics estimation during reception of information 
data sequences. 

21 . Method in accordance with claim 20 characterized in that for a data sequence, the interpolation is a linear interpo- 
lation performed starting from the channel characteristics estimation for a probe sequence immediately preceding 

45 and a probe sequence immediately following the data sequence. 

22. Apparatus for transmission and reception of digital signals comprising information data propagated in a communi- 
cation channel characterized in that it comprises a transmitting apparatus (1 1) and a receiving apparatus (14) at 
the two ends of the channel; with the transmitting apparatus in turn comprising the following: 

so 

an intercalating device (12) receiving at input the information data to be transmitted and outputing a frame 
made up of information data sequences to be transmitted intercalated with probe sequences consisting of 
known data; 

a transmitter (13) for sending the frame over the channel with the receiving apparatus comprising in turn: 
'■s a receiver (15) for reception of the frame from the channel; 

- a discriminator (1 6) separating the probe sequences and the data sequences from the received frame; 

- an estimator (17) receiving at input the probe sequences and supplying at output the channel characteristics 
estimation for the received probe sequences; 
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an interpolator (18) receiving at input the estimation output by the estimator and producing by interpolation the 
desired channel characteristics estimation for the data sequences; and 

an equalizer (19) receiving the estimation output by the interpolator and equalizing on the basis thereof the 
data sequence extracted from the discriminator to output the correct transmitted data sequence. 



JNSDOCIO: <EP 09i20l9A2l.> 



9 



EP 0 912 019 A2 



a) 



probe 


data 


probe 


data 


probe 


data 


probe 


data 


s 


D 


s 


D 


S 


D 


s 


D 



frame , block . t 



block 



b) 





probe 
Spre 


data 
D 


probe 
Spost 




probe 
Spre 


data 
D 


probe 
Spost 


1 


block 


i 




block 


i 



c) 



probe 


data 


probe 


data 


probe 


data 


probe 


Spre 


D 


Smid 


D 


Smid 


D 


Spost 



block | t 

, block , 



Fig.1 



3NSOCCID: <EP 0912019A2_I_> 



10 



EP 0 912 019 A2 




g*post' 







9 A 

..< 


^1 

^ i f 

pre->^ 


<-9 A 


post 

► 


probe 


data 


probe 


data 


probe 


t 




Ldata 1 Lpost 


i 



i_pre i_post 

Fig3 



NSCOCIO. <EP 091201QA2_I_> 



11 



EP 0 912 019 A2 




3NSDOCID: <£P 0912019A2_L> 



12 



« ri ii 11 ii in ill in, I iii 

(11) EP 0 912 019 A3 

(12) EUROPEAN PATENT APPLICATION 

(88) Date of publication A3: ( 51) lnt C , 7 : H04L 25/02, H04L 25/03 

02.10.2002 Bulletin 2002/40 

(43) Date of publication A2: 

28.04.1999 Bulletin 1999/17 

(21) Application number: 98203496.9 



(22) Date of filing: 19.10.1998 



(84) Designated Contracting States: 


(72) 


Inventor: Cusani, Roberto 


AT BE CH CY DE DK ES Fl FR GB GR IE IT LI LU 




00152 Roma (IT) 


MC NL PT SE 






Designated Extension States: 


(74) 


Representative: Faraggiana, Vrttorio, Dr. ing. 


AL LT LV MK RO SI 




(ngg. Guzzi & Ravizza S.r.l. 






Via Vinoenzo Monti 6 


(30) Priority: 22.10.1997 IT MI972376 




20123 Milano (IT) 


(71) Applicant: Telital S.p.A. 






34010 Sgonlco (Trieste) (IT) 







(54) Probe sequences intercalating method with information data for the estimation of 
communication channel characteristics and correct reception of data sequences 




Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



(57) A method for estimation of a transmission chan- 
nel and for transmission and reception of digital signals 
over she channel comprises the steps of transmission of 
data sequences intercalated with probe sequences, cal- 
culation m reception of a channel characteristics esti- 
mation lor probe sequences received and interpolation 
ol the channel estimation for the data sequences. After 
obtaining the interpolated estimation the transmission 
and reception method comprises the additional phase 
of recognition of received information data by equalizing 



the data sequences with the interpolation result. An ap- 
paratus employing the method has at one end of the 
chan nel a device (1 2) for intercalation of probe sequenc- 
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end a discriminator (1 6) which separates the probed and 
data sequences. An estimator (1 7) estimates the chan- 
nel for the probe sequences and an interpolator (1 8) pro- 
duces the desired interpolated estimation. An equalizer 
(19) equalizes the data sequence on the basis of the 
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